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Abstract We study the generation of a stochastic gravitational wave (GW) 
background produced by a population of neutron stars (NSs) which go over a hadron- 
quark phase transition in its inner shells. We obtain, for example, that the NS phase 
transition, in cold dark matter scenarios, could generate a stochastic GW background 
with a maximum amplitude of /ibg ~ 1CP 24 , in the frequency band f ^ s ~ 20—2000 Hz 
for stars forming at redshifts of up to z ~ 20. We study the possibility of detection 
of this isotropic GW background by correlating signals of a pair of Advanced LIGO 
observatories. 

Keywords Neutron Stars ■ Gravitational Waves • Strange Matter • Background 
1 Introduction 

The first scientific runs on the gravitational wave (GW) observatories have already 
been done and the first results have been obtained. In particular, we take the american 
interferometer LIGO that has already put some limits on the properties of neutron star 
(NS) ellipticities [T][2], for example. Putting some lower limits on such properties, like 
ellipticity, LIGO, among others like Virgo, GEO600 and TAMA, have contributed on 
the better understanding of NS properties. 

NSs are excellent candidates to the first detection of GWs, with black holes (BHs), 
due to its great amount of matter constrained into a very small radius. NSs are 
astrophysical objects on which a great amount of quadrupole moment variation can be 
generated by a mountain on the star, a crust quake, a core quake, a glitch, rotation or 
collapse. All these events would generate and radiate GWs. 

In this work, we address special attention to the GWs generated after a micro- 
collapse of the NS. Such event occurs, in a catastrophic way, when the NSs suffer a 
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phase transition in its interior. This phase transition, in this work, is assumed by the 
releasing of quarks from the interior of hadrons. When the critical density is reached, 
quarks deconfine. As the quark matter equation of state is softer then the hadronic 
one, the NS shrinks. 

It is widely accepted that NSs are born with fast rotating angular velocities. Due 
to magnetic torques, however, the NS periods could well be spun down. This spin- 
down causes a reduction in the centrifuge force and, consequently, the central energy 
density of the NSs increases. Those stars, born with densities close to that of the quark 
deconfinement, may undergo a phase transition forming a strange quark matter core. 
It is worth stressing, however, that it is not a common sense that such a transition 
really takes place. 

As a consequence of such a putative phase transition, the stars could suffer 
a collapse which could excite mechanical oscillations. The great amount of energy 
generated in this process, AE ~ 10 53 ergs (~ O.IMqc 2 ) [3J, could be dissipated, at 
least partially, in the form of GWs, which are studied in the present work. 

The energy released to excite these mechanical oscillations can be driven into many 
different modes, such as the fundamental f-mode, the pressure p-mode, which are the 
overtones of the fundamental mode, g-modes, excited by a sharp difference in the 
density distribution, the w-mode, which is a pure relativity mode. The energy is also 
driven to heat the star, releasing neutrinos or ejecting the crust material [4]. 

The number of NSs that undergo a phase transition was estimated to be about 
10 -6 events /year /galaxy [5]. However, this number depends on the properties of 
the phase transition and can increase (drop) to more (less) significant values. We 
investigate, in this work, the formation of a stochastic background of GWs originated 
by such event in scenarios of cosmological structure formation. 

We have improved the calculation made by Sigl[S] by using a realistic equation of 
state for describing NS matter, previously studied with politropic equations of state. 
For the rate of events, we have used the results obtained by Marranghello el. al. [5]- 
Finally, the amount of energy released in the GW mode is also estimated by realistic 
calculations made in ref. [7|. Moreover, we adopted in the present paper the history of 
star formation derived by [8], who employed hydrodynamic simulations of structure 
formation in a A cold dark matter (/1CDM) cosmology. A larger set of variables, such 
as the initial mass function and progenitor masses, is also used for a richer analysis. 
Even though the results do not show great differences, once the physics included is 
more realistic, the model becomes more reliable. 

Springel & Hernquist [8j obtain the history of star formation from hydrodynamic 
simulations of structure formation in ACDM cosmology. They study the history of 
cosmic star formation from the "dark ages", at redshift z ~ 20, to the present. They 
take into account besides gravity and ordinary hydrodynamics, radiative heating and 
cooling of gas, supernova feedbacks, and galactic winds. Their paper improves previous 
studies which consider either semi-analytical models or numerical simulations. 

It is worth mentioning that the story of star formation they obtain is consistent 
with observations. It is important to bear in mind, however, that nowadays observations 
give information from at most a redshift around z ~ 5. In the future, however, with 
the Next Generation Space Telescope (NGST) it will be possible to trace the cosmic 
star formation rate out to z > 20 (see, e.g., [9])- 

Besides the reliable history of star formation by Springel & Hernquist, we consider 
the role of the parameters 0:4=1,2, which gives the fraction of the progenitor mass which 
forms the remnant NSs. We consider that the remnant mass is given as a function of 
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the progenitor mass M r = a\ra + a.2- Recall that a given initial mass function (IMF) 
refers to the distribution function of the stellar progenitor mass, and also to the masses 
of the remnant compact objects left as a result of the stellar evolution. In section 3 we 
show in detail why M r is written as above. 

In the present study we have adopted a stellar generation with a Salpeter IMF, 
which is consistent with Springel & Hernquist, since they show that population II stars 
could have been formed at high redshift too. We then discuss what conclusions would 
be drawn whether (or not) the stochastic background studied here is detected by the 
forthcoming GW observatories such as LIGO and VIRGO. 

The paper is organized as follows. Section 2 deals with the NS equations of state 
adopted in our study; in Section 3 we consider the IMF and the NS masses; Section 
4 deals with the GW production; in Section 5 we present the numerical results and 
discussions; in Section 6 the detectability of the background of GWs are considered 
and finally in Section 7 we present the conclusions. 

2 Neutron stars and the nuclear matter equation of state 

The global properties of mass and radius of compact stars are directly related to 
its equation of state. In the present study we do not adopt a polytropic equation 
of state, since this could result in a poor modeling. In particular, we consider a model 
developed by |10j for a new class of parameterized field-theoretical model described 
by a Lagrangian density where the whole baryon octet is coupled to scalar and vector 
fields through parameterized coupling constants. 

The main characteristic of this model is to adjust these parameterized coupling 
constants to fit both the nuclear matter properties and the neutron stars observation 
data. Varying this parametrization, we are also able to fit the properties of some other 
field-theoretical models like the Walecka model and the Zimanyi-Moskowski model and 
compare our results to those of other relevant models in the literature as described by 
[10] . We are also able to reproduce calculations with higher and lower values for neutron 
star masses and radii, compression modulus and effective baryon masses. 

Instead of performing the following calculations with a wide range of equations 
of state we preferred to describe our work choosing a representative equation of 
state whose coupling constants were chosen to be those that fit both the nuclear 
and stellar properties in an acceptable manner. It is important to notice that other 
parameterizations of this model could result in quite similar results, harder detection 
of the background or unreliable results due to the not well fitted characteristics of 
nuclear matter. 

The nuclear matter Lagrangian density is described by 

c = ^2 ^ B (^(^ ~~ oZb^) - ( m b - g^B a ) - [\g* B B~r ■ el) 1>b 

B 

i / o nft 2 2\ 1 uv . 1 2 n 

+ -(a M (Td f a — m a a ) — -lu^lu^ + -itl^uj^lu^ 

-\q^ ■ + \ m2 gQ^ ■ Q 1 " 1 + ^2 4>i[ilnd^ - Mi]ipi , (1) 

/ 

where the baryon field ipg is summed over the whole baryon octet and coupled to the 
scalar and vector fields a, lu and g through the parameterized coupling constants gig, 



Inn 



300 400 
K (MeV) 



Fig. 1 The NS maximum masses versus the compression modulus, for different 
parametrization in the Taurines et.al. models. The solid line represents a parametrization 
in the scalar mesons while the dotted line describes a parametrization in both scalar and 
vector contributions of the nuclear force. 



<?* B , <7gg. The free lepton fields ip\ contributes to the electrical equilibrium in the NS 
matter. The masses of the baryons, mesons and leptons are represented by Mg, m CT ,ct>,@ 
and m\, respectively. The free lepton fields contributes to the electrical equilibrium in 
the NS matter. 

Using a parametrization for the baryon-meson coupling constants, 
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this model describes a wide range of NS parameters, such as a maximum mass ranging 
from (very low values) M — 0.66Mq up to M = 2.77 Mq, and the corresponding radii 
that vary in the range of 8 < R < 13 km. Each pair in the mass-radius relation is 
associated with a different parametrization of the equation of state (see [10] for further 
details). Obviously, each value of such parametrization represents different values of 
nuclear matter properties. 

The study of nuclear matter properties, however, is still a very open problem. A 
relevant parameter is the compression modulus, which defines the curvature of the 
equation of state or, in a few words, is related to the capability of matter to be 
compressed. One can see in Figure \T\ that the maximum mass of NSs is obtained 
for large values of the compression modulus. The definition of K reads 



K = 9 



2d 2 (e/p) 



dp 2 



(3) 



where e is the energy density, p is the baryon density and po the baryon saturation 
density. Actually, the acceptable values for the compression modulus are constrained 
to be in range 200 < K < 300MeV [TTHT3] . 

An important constraint to the equation of state is to describe a family of NSs 
with a maximum mass greater than the value obtained by the most massive known 
pulsar. The PSR J0751+1807 has a mass M = 2.1 ± 0.2M Q [13]. This condition has 
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to be fulfilled by the set of parameters chosen to describe the NS model. This first 
condition already constrains the value of compression modulus to K > 22QMeV when 
this model is considered. Smaller values of K would lead to masses smaller than 1.9Mq 
disagreeing with the measures of PSR J0751+1807 mass. 

It is worth mentioning that quite recently [TJ reported the discovery of eight new 
millisecond pulsars in NGC 6440 and NGC 6441. In particular, they could determined 
the mass of the pulsar J1748-2021B (NGC 6440B), namely M = 2.74±0.21M Q . If the 
mass of this pulsar is really high, this would mean that the value of the compression 
modulus would be around 500MeV . In other words, our model could describe, in an 
extreme limit of coupling constants, hybrid stars with that value of gravitational mass. 
However, it would lead to values of nuclear matter properties that are not in complete 
agreement with experiments. 

Another scenario to be discussed concerns the strange quark star. The model used 
in this analysis is the MIT bag model [15] • The MIT bag model describes, as its own 
name suggests, that baryons are composed by three quarks confined inside a bag. 
The effects of pressure difference in the interior and exterior regions of the bag are 
summarized in the bag constant. 

Even with considerable advances in the study of QCD (quantum chromo-dynamics) 
on the lattice, from where we expect the most reliable theoretical results and the 
consequent advances in collision experiments, the bag constant still presents a very 
wide range of possible values. In the same way, the perturbative constant, a c , the 
one that represents the first order correction to the strong interaction forces between 
quarks, also represents an open issue. 

Once the transition to quark-gluon matter occurs, the weak interaction processes 
for the quarks u, d and s 

u + s — > d + u (4) 

and 

d + u — > u + s (5) 

will take place, and a rapid transition occurs with a consequent gravitational micro- 
collapse. As a result, a huge amount of energy is dissipated, some in the form of GWs 
produced by quasi-normal modes excitation. The amount of energy released during the 
transition is calculated to be around 10 53 erg according to ref. [3]. Only a small part of 
this energy is released in the form of GWs [7]. 



3 Initial mass function and the neutron star masses 

The calculation of the GW background from NS phase transition requires the knowledge 
of the distribution function of stellar masses, the so called stellar initial mass function 
(IMF), (f>(m). Here the Salpeter IMF is adopted, namely 

4>(m) = Am- (1+X \ (6) 

where A is the normalization constant and x = 1.30 (our fiducial value). The 
normalization of the IMF is obtained through the relation 



/>m u 

j m<j>(m)dm = 1, 

mi 



(7) 
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Table 1 The values of 0^=1,2 which determines the fraction of the NS remant. 



~a x 1732 I75O 17T7 

a 2 (M ) 3/4 1.1 0.53 



where we consider m\ = 0.1 Mq and m u = 125 Mq. For further details we refer the 
reader to, e.g., [16| . 

It is worth mentioning that concerning the star formation at high redshift, the IMF 
could be biased toward high-mass stars, when compared to the solar neighborhood IMF, 
as a result of the absence of metals [171118] . 

On the other hand, in the study by Springel and Hernquist, besides a high- mass 
star formation, it is shown that the population II stars, whose IMF could well be of 
the Salpeter's type, could start forming around redshift 20 or higher. 

In the present study we consider this population II studied by these authors. Then, 
for the standard IMF, the mass fraction of NSs produced as remnants of the stellar 
evolution is 



where m m i n is the minimum stellar mass capable of producing a NS at the end of its 
life, and M r is the mass of the remnant NS. Stellar evolution calculations show that the 
minimal progenitor mass to form NSs is m m j n = 8Mq , while the maximum progenitor 
mass is m ma x = 25 — 40Mq (see, e.g., [H]). 

We do not intend to discuss in the present paper stellar evolution scenarios related 
to the formation of NSs. Many works can be found in the literature concerning the study 
of stellar evolution, supernova and NS formation 20 . However, we need to formulate 
a relation between the progenitor star mass m and the remnant NS mass M x . For the 
remnant, M r , we take 



where ol\ and 02 are constants; the first one is dimensionless and the second one is 
given in solar masses. As the results for stellar evolution are not yet fully determined, 
we studied this parametrization in three different scenarios for the values of ai and a 2 
which would represent the distribution of NS mass as function of the progenitor mass 
(see Table 1 and the following sections). For example, considering NSs formed from 
progenitors with m — 20Mq we would have remnants with masses M r = 1.375M0, 
1.5M© and 2.7 Mq. These parameters describe sharper or softer distribution of NS 
masses around the standard value of mass M = IAMq . Note that we are considering 
in the end that M r is independent of the redshift. 

With these considerations at hand, the mass fraction of NSs reads up to /jjg = 10~ 2 
for x — 1.30, while the fraction of NSs that undergoes a phase transition (/|jg) can 
drop down to values <C 1 (see, e.g., [5]). 

To assess the role of possible IMF variations in our results, other values of x have 
also been considered. Besides the standard IMF, two others have been studied, namely, 
with x — 0.3 and x = 1.85, which yield approximately ten times and one-tenth of mass 
fraction of NSs of the standard IMF, respectively. 




(8) 



M r = Qi m + Q2 , 
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Table 2 Models description: mass function constants, ai,2, progenitor, m, and remnant, M, 
mass ranges, redshift, z, and the observed frequency, v. 



Model 


<X\ 


Q-2 


Am(M Q ) 


AM(Mq) 


Az 


Au(Hz) 


A 


1/32 


3/4 


8.0-40.00 


1.00-2.00 


0-20 


48-2000 


13 


1/32 


3/4 


14.4-16.32 


1.20-1.26 


0-50 


31-1666 


c 


1/32 


3/4 


14.4-16.32 


1.20-1.26 


0-20 


75-1666 


D 


1/32 


3/4 


20.8-22.72 


1.40-1.46 


0-20 


65-1428 


E 


1/50 


1.1 


8.0-40.00 


1.26-1.90 


0-20 


50-1587 


F 


1/50 


1.1 


15.0-18.00 


1.40-1.46 


0-20 


65-1428 


G 


1/17 


0.53 


8.0-25.00 


1.00-2.00 


0-20 


48-2000 



In order to cover a wide number of parameters we present, in Table 2, the models 
considered in our study. In the first column appears the name of the model; in the 
second and third columns we present different combinations of the parameters ot\ and 
lv2j which imply in different ways to calculate the NS remnant mass for a given IMF; 
in the fourth column the mass range of the progenitor star; in the fifth column the 
NS remnant mass; in the sixth column the NS redshift formation, and finally in the 
seventh column the observed frequency, which are obtained via the empirical formulae 
given in reference |21] , 

These parameters are directly related to the number of NSs that goes over the 
phase transition. According to the NS models obtained by realistic equations of state, 
only a fraction of them develops a core composed by deconfined quarks and gluons. 
When we set the mass range of the progenitor mass, we are also setting the mass 
of the remnant object. Stars with higher mass reach central densities that are high 
enough to develop the quark core, while smaller stars do not. The available energy 
for GW emission generated in the collapse process is also implicit in the mass of the 
remnant star. The more massive is the star, greater is its core and greater is the energy 
difference between the neutron and the hybrid star. We have found that stars with 
baryonic mass M(, = 2Mq develops a core as large as R — 7km and can generate as 
much as 3 x 10 53 erg of energy, while stars with baryonic mass Mf, = IMq develops a 
small core with R = lkm and has ten times less energy available from the transition. 
As we do not know for sure the amount of this energy that is driven in each mode and 
in order to simplify our calculations we have adopted the standard medium value of 
O.OIMq for the released energy. 



4 Gravitational wave production 

The GWs can be characterized by their dimensionless amplitude, h, and frequency, 
v. The spectral energy density, the flux of GWs, received on Earth, F v , in 
erg cnT~ 2 s -1 Hz -1 , is (see, e.g., [221125] ). 

where lo q -^ s = 27ri/ obs , with the GW frequency (Hz) observed on Earth, c is the 
velocity of light, G is the gravitational constant and is the strain amplitude of the 
GW in Hz -1 / 2 . 
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The stochastic GW background produced by NS phase transition would have a 
spectral density of the flux of GWs and strain amplitude also related to the above 
equation (|9]). Therefore, in the above equation the strain amplitude takes into account 
the star formation history occurring at the 'first light', just after the 'dark age' epoch. 
The strain amplitude at a given frequency, at the present time, is a contribution of NSs 
with different masses at different redshifts. Thus, the ensemble of NSs formed produces 
a background whose characteristic strain amplitude at the present time is y^h- 

On the other hand, the spectral density of the flux can be written as [24,25 

rz c \ rm u 

Fu= / M^obs^NsK*), (io) 

J z c { J m m i n 

where /i/(f bs) i s the energy flux per unit of frequency (in erg cm _2 Hz _1 ) produced 
by the formation of a unique NS and di?NS i s the differential rate of NSs formation. 

The above equation takes into account the contribution of different masses that 
collapse to form NSs occurring between redshifts z c ; and z c { (beginning and end of the 
star formation phase, respectively) that produce a signal at the same frequency v h s . 
On the other hand, we can write /i/(^ bs) EB] as 

7TC 3 2 

M^obs) = "^T ft NS. (11) 

where /j-ns is the dimensionless amplitude produced by the NS micro-collapse of a given 
star with mass m that generates at the present time a signal with frequency ^obs- Then, 
the resulting equation for the spectral density of the flux is 

™ 3 f u2 



F »=2G I h z NS dR NS . (12) 



From the above equations we obtain for the strain amplitude 



1 



s h = / /4 s di? NS . (13) 

^obs 



Thus, the dimensionless amplitude reads 

h% G = — — / /in S gLR ns , (14) 

^obs J 

(see [271128]). 

The micro-collapse produces GWs at frequency v of the NS f-mode, and 
dimensionless amplitude given by |29| 

, ,„-19 f E \ 1/2 (2kHz\ 1 / 2 flMpc\ 

^ lxl ° \m^) {—) l^t) (15) 

where E is the available pulsation energy, and is the luminosity distance to the 
source. 

Recall that the energy available to excite the pulsating modes is directly related to 
the equation of state adopted. 

Note that, since we are considering cosmological sources we have to take into 
account the redshift effect on the emission frequency, that is, a signal emitted at 
frequency v c at redshift z is observed at frequency i/ b s = "e(l + z)^ 1 . 
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For the differential rate of NS we consider only those that undergo micro- collapse, 
namely 

dRm = p*(z)-^-f^ s cf>(m)dmdz, (16) 

where p*(z) is the star formation rate (SFR) density (in M^yr -1 Mpc -3 ), dV is the 
comoving volume element, and /j^g, as already mentioned, is the fraction of NSs which 
may undergo a phase transition forming a strange quark matter core. 
The comoving volume element is given by 

dV = iir(jj-y*F(n M7 n A ,z)dz, (17) 



with 



^(D M ,D A ,z) = 7=== == = =, (18) 

+ ^) 2 (i + n M z) - z{2 + z )n A 



and the comoving distance, r z , is 




where 



Jo F{D M ,D A ,z') J 



Dm = fl DM + and 1 = D k + + D A (20) 



are the usual density parameters for the matter (M), i.e., dark matter (DM) plus 
baryonic matter (B), curvature (k) and cosmological constant (A). The function S is 
given by 

{sin x if closed, 
x if flat, (21) 

sinhx if open. 

The comoving distance is related to the luminosity distance by 

d L =r z {l + z). (22) 

The set of equations presented above can be used to find the dimensionless 
amplitude of the GW background produced from a population of NSs which goes over 
a hadron-quark phase transition in its inner shells as a function of the SFR density, 
and related to the 'first light' epoch. 

For the SFR density, we adopt the one derived by Springel and Hernquist (see, [8] 
for details), namely 

. , , _ f3e*p[A{z- z m )] , . 

(3 - A + A exp [f3(z - z m )\ 

where A = 3/5, /3 = 14/15, Zm — 5.4 marks a break redshift, and pm — 
0.15 MQyr _1 Mpc -3 fixes the overall normalization. 

It is worth mentioning that these authors employed hydrodynamic simulations of 
structure formation in a TlCDM cosmology with the following parameters: Dyi = 0.3, 
D A = 0.7, Hubble constant H = 100 h kms -1 Mpc -1 with h = 0.7, D B = 0.04, 
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and a scale-invariant primordial power spectrum with index n = 1, normalized to the 
abundance of rich galaxy clusters at present day (erg = 0.9). 

Another relevant physical quantity associated with the GW background is the 
closure energy density per logarithmic frequency span, which is given by 

%W = -^. (24) 
The above equation can be rewritten as 

%w = = ^l hs hl G . (25) 

In the next section we present the numerical results and discussions, which come 
mainly from the equation for /ibg- 



5 NUMERICAL RESULTS AND DISCUSSIONS 

In this section we present the numerical results of the putative GW background related 
to the NS phase transition. 

The Eq. 15 is the main equation to calculate such a GW background, with the 
relevant ingredients given by Eqs. 16 — 24. Basically, to integrate Eq. 15 one needs 
to specify the cosmological model, the star formation rate density, etc. We refer the 
reader to Section 4 to see in detail how the various relevant ingredients are taken into 
account. 

More explicitly, to obtain the spectrum of the background using Eq. 15 one takes 
a given v Q h s and perform integrations in the variables m and z, where the integration 
intervals are those of the models present in Table 2. This procedure is repeated for 
different f bs' s i whose values are taken from the frequency band also present in Table 
2. 

In Figure [2] we compare the background of GWs generated by NSs, which undergo 
phase transition, with that by the BH formation, which undergo quasi-normal mode 
instability (see [30] for details). In particular, for the NSs we consider two situations. 

We consider in Figure 2 two versions of the model B (see Table 2): one with its 
original values, namely, /^g ~ 0.01 and the other one, just as an example, with /™ = 1. 
Obviously, even the most favorable values for nuclear and subnuclear matter coupling 
constants could not make even tend to 1, but this model is included just as a 
matter of comparison, as it would represent an upper limit. In both cases we consider 
E ~ O.OIMqc 2 , which means that approximately only a tenth of the energy generated 
in phase transition goes to the pulsating mode. Note that we are considering an energy 
release in GW about ten times higher then [6]. Later on we comment on how this 
choice modifies our conclusions concerning the putative detectability of the background 
studied here. 

Note that the background of GWs generated by the NSs would have an amplitude 
greater than that generated by the BHs only if a considerable fraction of the NSs 
formed undergo phase transition. 

In Figure |3] we compare how different laws to calculate the NSs mass, the remnant 
mass, for a Salpeter IMF, modify the spectrum of the background of GWs generated. 
Note that there is no significantly difference in the background for the three cases 
considered. 
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Fig. 2 The dimensionless amplitude for the GW background of stellar BH formation (solid 
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Fig. 3 The dimensionless amplitude for the GW background of micro-collapsed NSs for 
1.2 < M < 1.26M (solid line), 1.4 < M < 1.46M (dotted line) and a sharper NS 
distribution with 1.4 < M < 1.46Mq (dashed line). We consider in these models /^* g = 0.01 
and E = 0.01M Q c 2 
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Fig. 4 The dimensionless amplitude for the GW background of micro-collapsed NSs for a 
IMF with x = 1.30 (solid line), 0.30 (dotted line) and 1.85 (dashed line). We consider in these 
models f*L = 0.01 and E = O.OlM c 2 . 



We also investigated the role of the IMF in our results. In particular, besides the 
standard IMF (x = 1.30) we also consider two others, namely, x — 0.30 and 1.85. This 
same choice was considered by [30] in their study on the background of GWs generated 
by the formation of population III stellar BHs. As with for the BHs the amount of NSs 
increases (decreases) for x = 0.30 (1.85) with respect to x — 1.30. In Figure |4] we 
compare these three cases. In this comparison, we consider model B. 

We now compare model B, where we consider ~ 0.01 and E ~ 0.01 Af Q c 2 , 
with the spectra of backgrounds of GWs for two other pulsating modes, considering 
these two very parameters. In particular, we present in Figure the GW background 
of micro-collapsed NSs considering the contributions of the p-mode at 6 kHz (dotted 
line) and the w-mode at 12 kHz (dashed line). As expected, see Figure(5] the higher the 
frequency is, the more the spectrum is shifted to the higher frequencies. Concerning 
the amplitude, the spectra do not present considerable differences among them. 

It is worth mentioning that we did not take into account the g-mode in our 
calculations because the amount of energy put into this mode, which is related to 
density discontinuities [4], is much smaller then in the case of f- and p-modes and the 
frequencies related to this mode lie in a lower range of the spectrum. 

A relevant question is whether the background we study here is continuous or not. 
The duty cycle indicates if the collective effect of the bursts of GWs generated during 
the collapse of a progenitor star generates a continuous background. The duty cycle is 
denned as follows: 




(26) 
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Fig. 5 The dimensionless amplitude for the GW background of micro-collapsed NSs 
considering the contributions of different quasi- normal modes. The f-mode (solid line), the 
p-mode at 6kHz (dotted line) and the w-mode at 12kHz (dashed line). We adopted model B 
(f-mode) with = 0.01 and E = O.OIMqc 2 . For the p- and w- modes we also considered for 
comparison these very two parameters. 



where Atq^j is the average time duration of single bursts at the emission, which 
amounts to ~ 100 ms to the f-mode of NSs, as obtained though the empirical equations 
of ref. |21J . This value for the f-mode damping time is quite different from that adopted 
by Sigl [BJ, who assumed a value of 1/3 ms. 

Since the star formation rate could be high, a significant amount of GWs could be 
produced. We also note that, independently of the primordial cloud mass and of the 
redshift of collapse, star formation occurring at high redshift could produce high duty 
cycle values, which lead us to conclude that the stochastic GW background could be 
continuous. For all the models studied here the duty cycle is > 1. 

In the next section we consider the detectability of the background of the GW 
background that we propose exists. 



6 DETECTABILITY OF THE BACKGROUND OF GRAVITATIONAL 
WAVES 

The background predicted in the present study cannot be detected by single 
interferometric detectors, such as VIRGO and LIGO (even by the forthcoming 
advanced ones). However, it is possible to correlate the signal of two or more detectors 
to detect the background that we propose exists. Christensen [31] was the first to show 
that this kind of signal can, in principle, be detected by correlating the outputs of two 
different detectors. However, the main requirement that must be fulfilled is that they 
must have independent noise. This study was improved by |32| and by [33]. The reader 
should also refer to the papers by [34] and |35] who also deal in detail with such an 
issue. 
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Table 3 For the models of Table 2 with different values of /jy S , we present the S/N for pairs of 
Initial, Enhanced and Advanced LIGO observatories for one year of observation. We consider 
in the calculations E = O.OIMqc 2 . 





■Ins 




S/N 




Model 


Initial LIGO 


Enhanced LIGO 


Advanced LIGO 


A 


1.0 


8.6 x 10" 4 


2.0 x 10~-> 


2.4 x 10 _i 


B 


0.01 


2.2 x 10~ 6 


4.8 x 10~ 6 


2.0 x 10~ 4 


C 


0.01 


2.2 x 10~ 6 


4.8 x 10~ 6 


2.0 x 10~ 4 


D 


0.01 


2.8 x 10" 6 


6.3 x 10~ 6 


2.6 x 10" 4 


E 


1.0 


1.1 x 10" 4 


1.5 x 10~ 3 


3.0 x 10~ 2 


F 


1.0 


5.6 x 10~ 6 


1.0 x lO" 5 


5.7 x 10~ 4 


G 


1.0 


6.5 x 10" 4 


1.5 x 10" 3 


7.6 x 10~ 2 



To assess the detectability of a GW signal, one must evaluate the signal-to-noise 
ratio (S/N), which for a pair of interferometers is given by (see, e.g., [©,34 ) 

(S/N) 2 = 

where S^' is the spectral noise density, T is the integration time and is the overlap 
reduction function, which depends on the relative positions and orientations of the two 
interferometers. For the j(v) function we refer the reader to [33], who was the first 
to calculate a closed form for the LIGO observatories. Flanagan [33] (see also [34]) 
showed that the best window for detecting a signal is < v < 64 Hz, where the overlap 
reduction function has the greatest magnitude. 

Here we consider, in particular, the LIGO interferometers: the Initial, the Enhanced 
and the Advanced ones. For the Initial LIGO, the present configuration, the sensitivity 
curve has been taken from |http://www.ligo.caltech.edu/^jzw eizig/distributi on/| 
LSCJData/. It is worth mentioning that the design-sensitivity curve was almost 
attained by the S5 run. For the Enhanced LIGO, which is an upgrade of 
the Initial LIGO and to be operative briefly, the sensitivity curve has been 
taken from http://www.ligo.caltech.edu /docs/T/T060156-01.pdf, in particular from 
its Fig. 12. The Enhanced LIGO is approximately a factor of two more 
sensitive than the Initial LIGO. Finally, for the Advanced LIGO, designed 
to be approximately a factor of ten more sensitive than the Initial LIGO, 
the sensitivity curve has been taken from http://lhocds.ligo- wa.caltech.edu:8000 
/advligo/LSC_Modeling?aetion=AttaehFile&do=get&target=isccddFebl9.pdf, in 
particular from its Fig. 1. 

We show in Table 3 the S/N for the models of Table 2 for the three different LIGO 
generations and for different values of ffjg- 

As shown in Table 3, the signal-to-noise ratio for all models studied is lower than 
one, even for an advanced LIGO. Therefore, contrary to the claim of [6] such a putative 
GW background would hardly be detected. 

Using the same parameters of model B with = 0.01 and E — O.O1M0C 2 , we 
have also investigated the role of the variations in the IMF modifying the exponent 
x. We have found, for the Advanced LIGO, a signal-to-noise ratio of 3.4 x 10~ 5 , for 
x = 0.30, and 1.5 x 10 -4 , for x — 1.85. As before, the signal-to-noise ratios are lower 
than 1. 



50tt 4 



T 



f/;/ 



(27) 
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We still have to mention the results for different frequencies, which correspond to 
the excitation of p- and w-modes. From our calculations, with the same parameters 
shown for model B , i.e., with /^g = 0.01 and E — O.OIMqc 2 , we have changed the 
frequencies to v e — 3, 6 and 12 kHz, we obtain a signal-to-noise ratio for Advanced 
LIGO equal to 2.8 x 10 -5 , 7.8 x 10 -6 and 1.3 x 10" 8 , respectively. These results show 
that such background would not either be detected. 

Note that the signal-to- noise ratio, for given IMF, aj, and integration time, depends 
on /j^g and E as follows 

(S/N)oc/^ s S; (28) 

and it also depends on the SFR density in a more complicated way, namely, through 
an integral involving the redshift z. The higher the star formation rate, the higher the 
signal-to noise ratio will be. 

Just as a matter of comparison, even considering models with /™ = 1, a signal- 
to-noise ratio significantly greater than one for the Advanced LIGO would be possible 
either the SFR density would be much greater than that by Springel & Hernquist, 
or the energy generated in the phase transition were almost completely channeled to 
excite the f- mode. Obviously, an optimist combination of the SFR density and the 
energy channeled to the f-mode would also render the same. 

Still concerning the energy release in GWs, the above equation shows explicitly 
how the detectability is affected for lower values of E. For example, adopting an energy 
release as that by [B] the signal-to-noise ratios presented in Table 3 would be ten times 
lower. 

7 CONCLUSIONS 

We present here a study concerning the generation of GWs produced from a 
cosmological population of NSs. These stars may undergo a phase transition if born 
close to the transition density, suffering a micro-collapse and exciting quasi-normal 
modes. 

This work differs from the previous one done by Sigl in two essential points. First, 
we have calculated the NS star structure using realistic equation of state, instead 
of a polytropic model, obtained through a field-theoretical model which takes into 
account the microstructure of matter composed by the whole baryon octet, coupled 
to meson fields. We have also taken into account, into the quark sector, the standard 
MIT bag model. The transition is constructed via Gibbs criteria. We have calculated 
the amount of energy released by the transition and the rate of stars that undergo this 
microcollapse. Second, we have used another scenario for the cosmical star formation 
rate. We also have made a study considering various scenarios varying the set of 
parameters that gives the main characteristics of star formation masses, etc. 

We show that a detectable background is possible only if the SFR density is much 
greater than that predicted by Springel & Hernquist or if the energy generated in the 
phase transition is almost completely channeled to excite the f- mode. Obviously, a too 
optimistic combination of these two possibilities could do the same. 

It is worth mentioning that even planned detectors like DECIGO and BBO would 
not be sufficiently sensitive to detect the background considered here. 

Finally, it is worth mentioning that Sigl claims a putative marginal detection if 
the majority of NSs undergo a phase transition. In our study we obtain a similar 
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conclusion for some of our models if the above condition is fulfilled (/jjg = 1) and for 
an integration time significantly larger than one year. 

One could ask why we obtained such a conclusion even adopting a release of GW 
energy larger than that adopted by Sigl. One reason for that is related to the SFR 
densities we adopted, namely, in his case a larger number of NSs is generated. As a 
matter of comparison, if you multiplied the SFR density by a factor of ten the signal- 
to-noise ratio would be multiplied by this very factor. Other reasons could be related 
to the way we calculate the NSs remnants (see Table 3) and the value Sigl adopted for 
the damping time. A damping time of ~ 100 ms instead of one of 1/3 ms could modify 
his results significantly. 
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